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We compute the universal generic corrections to the power spectrum in slow-roll inflation due to 
unknown high-energy physics. We arrive at this result via a careful integrating out of massive 
fields in the "in-in" formalism yielding a consistent and predictive low-energy effective description 
in time-dependent backgrounds. The density power spectrum is universally modified at leading 
order in H/M, the ratio of the scale of inflation to the scale of new physics; the tensor power 
spectrum receives only subleading corrections. In doing so, we show how to make sense of a physical 
momentum-cut-off in loop integrals despite dynamical redshifts, and how the result can be captured 
in a combined effective action/effective density matrix, where the latter contains non-adiabatic terms 
which modify the boundary conditions. 
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INTRODUCTION 

The dawning of the era of precision cosmology de- 
mands that we understand the history of our universe 
theoretically with the same accuracy as experiment. The 
WMAP determination of acoustic peaks in the CMB 
spectrum to 1% accuracy has given strong support to the 
existence of an era of inflation [T]. Inflation, famously, 
predicts the primordial power spectrum of density fluc- 
tuations underlying all structure in the universe. 

This ability of WMAP and future Planck data to 
constrain theoretical models of the early Universe has set 
off a scramble to delineate a theoretically controlled com- 
putation of the primordial inflationary power spectrum. 
The textbook approach makes a number of explicit and 
implicit assumptions, each of which can affect the power 
spectrum at the accuracy measured. The main assump- 
tion is that the density spectrum can be reduced to that 
of an adiabatic fluid whose excitations are always weakly 
coupled to gravity. This can be conveniently encoded in 
a single scalar field Lagrangian with Bunch-Davies initial 
conditions in a fixed inflationary cosmology without met- 
ric fluctuations. The most obvious issue with these as- 
sumptions is that the redshifts implied by the 60 e-folds of 
inflation necessary to solve the horizon and flatness prob- 
lem, place the relevant momentum scales in an energy 
regime far beyond the Planck-scale, where in principle 
gravitational backreaction and quantum-gravity correc- 
tions cannot be ignored ;3L Turning the issue around, it 
also implies a window of opportunity that quantum grav- 
ity or any other New Physics arising at high energy scales 
could have a measurable effect on the power spectrum. 
This question was actively pursued some time ago with 
the conclusion that in toy models [3H52] one can obtain 
measurable corrections of the order H/M , comparable to 
intrinsic cosmic variance, with H the Hubble scale and 
M the scale of New Physics [55P0] . 



To truly connect with the data one needs the univer- 
sal generic model-independent corrections to the power 
spectrum in terms of an effective field theory that encodes 
order by order the corrections due to New Physics at high 
scales as well as deviations from adiabaticity. How to 
account for New Physics is formally well-understood in 
terms of Wilsonian effective actions, and for the adiabatic 
mode this has been actively pursued recently [2HJ PflT - 
146] [57] . However, Wilsonian effective actions are only 
consistent provided adiabaticity is maintained [47U49) . 
In practice this has meant that energy is assumed to be 
a conserved quantity. But precisely this is impossible in a 
cosmological time-dependent background. Redshifts con- 
tinuously mix the regimes of various scales and strictly 
speaking a well-defined separation of energy cannot be 
maintained. This long-standing paradox has fundamen- 
tally hampered the construction of low energy effective 
theories in cosmological spacetimes, literally since energy 
is not a conserved quantity. 

In a previous letter [50] we provided an algorithmic 
solution to this obstacle to compute the generic new 
physics corrections to the inflationary power spectrum. 
One can generate the universal low energy effective ac- 
tion by integrating out a massive field in any particular 
New Physics model. This is sensible in a cosmological 
setting, as long as one computes late time expectation 
values directly in non-equilibrium real-time QFT via the 
Schwinger-Keldysh approach. One of the principal dif- 
ficulties found in previous approaches was the correct 
procedure to implement an energy and momentum scale 
cutoff in a time-dependent background; in intermediate 
steps the order of energy and momentum integrals now 
no longer commute. We have overcome this through the 
realization that the contributions from the massive field 
may be reliably captured in a stationary-phase approxi- 
mation for vertex evaluation. This localizes interactions 
to fixed moments in time and allows us to unambigously 
implement the cutoff in momenta. 
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There are several noteworthy features that this algo- 
rithmic solution reveals: 



SINGLE FIELD SLOW-ROLL INFLATION PLUS 
A MASSIVE FIELD 



• Most importantly, New Physics and non-adiabatic 
corrections no longer separate. Integrating out the 
heavy field not only causes changes in the action, 
but also in the effective density matrix of the low- 
cncrgy field. The former can account for non- 
adiabatic effects. 



The non-adiabatic terms arc non-local in position 
space, but localized on a New Physics Hypersur- 
face where the physical momentum p(t) = k/a(t) 
equals the mass of the heavy field p(t) = M. In- 
tuitively this is what should happen and was the 
basis of many earlier ad hoc models. Here it is a 
consequence of integrating out the heavy field. 



• These terms can be interpreted as modified ini- 
tial conditions for the low-energy adiabatic infla- 
ton. This again confirms the qualitative insights 
gained from toy models, but now quantitatively. 

In this article we apply this cosmological effective field 
theory approach to the cosmologically relevant scenario 
of slow-roll inflation. We compute the corrections to the 
scalar power spectrum resulting from interactions with 
a heavy field. Parametrizing the corrections in a man- 
ner expedient to comparison with observational data, the 
universal generic correction to the scalar power spectrum 
due to unknown high-energy physics equals 



H 

M 



(Co + + 



where e\ is the first slow-roll coefficient. The power of 
the approach, however, is that each of the variables Cj 
can be computed in terms of the parameters of whichever 
theoretical model for the unknown New Physics one has 
in mind. For the power spectrum the result is merely an 
amplitude correction, but it will be very interesting to 
see the effect in higher-order correlation functions. 

This article is organized as follows. In §2 we present a 
simple model of slow-roll inflation containing new physics 
at high energies and analyze the field fluctuations in two 
gauges which will prove useful. In §3 we review the in-in 
formalism and how to compute field fluctuation correla- 
tions for our theory. In §4 we calculate the scalar and ten- 
sor power spectra in the spatially flat gauge, then convert 
these to the uniform density gauge and parametrize the 
answer in a way which facilitates easy comparison against 
observation. In §5 we discuss observational prospects and 
conclude. 



The minimal action description for inflation is a scalar 
field coupled to gravity, 



Sinful = / d 4 x^g 



M 2 pl R--(dcj > ) 2 -V {ni {4>) 



(1) 



We assume a spatially homogeneous ansatz for the back- 
ground metric 

ds 2 = -dt 2 + a(t) 2 dx 2 . 

Defining H = a/ a, the background equations of motion 
are 



12tt 



4-7T 

32^ 2 



v{4>). 



(2) 
(3) 



Primes denote derivatives with respect to the field, while 
overdots denote derivatives with respect to coordinate 
time. 

To facilitate analysis of inflating solutions, one defines 
the (first) Hubble slow-roll parameter ei 



M, 



pi 

47T 



H'{<t>) 
H(cf>) 



(4) 



Inflation occurs for t\ < 1 and ends when e\ — 1, which 
we use to define t on( j. It will prove useful to also define a 
second slow-roll parameter, 



62 fa) = 



dlnei 
dN 



M 2 \H"(<j)) 



One could continue defining an infinite hierarchy of such 
parameters e„, but this will be sufficient for our purposes. 

Of course there is a direct relationship between Vi n f(<j>) 
and ei, €2- This is given by 



V(4>) 

dV 
# 
d 2 V 
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= -3H 2 M n 



= ZH 
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It will be helpful to define the so-called "potential slow- 
roll" parameters 



ev(<P) 



M 2 , 
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i}v{4>) 
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We can then invert this to express ei,e 2 in terms of 
ey , rjv , with exact solutions possible in the case of power- 
law inflation. 
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Perturbative Interactions 

We may take the initial value of the inflaton to be 
4>(tin) = 0, so that the potential near this point is ap- 
proximated as 

where Vo,j, m 2 ,go are all constants. Thus ey,ryy, and 
hence e±, £2, can be solved for in terms of these coupling 
constants. 

To the inflationary action ([I]) we then add a massive 
field x with renormalizable interactions to the inflaton: 



S n 



■x\ 



d 4 x^/g 



1 

+ 2 
■0 2 X 
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4>x* 



We will be working to cubic order in field fluctuations, 
but are only interested in single contractions of the heavy 
field fluctuations, justifying our expansion in \ to second 
order. The combined action S = Si n { + S ncw will then 
produce an inflationary period but contain New Physics 
at scales determined by the couplings gi and mass scale 
M. 

A naive approach, based on static background quan- 
tum field theory intuition, might be the following: since 
we are presumed to be below the energy scale where \~ 
quanta can be created, we can simply integrate out this 
field to yield the following effective 'New Physics' poten- 
tial for d>: 



Vnew(<t>) 



1 



□ - M 2 



AM 2 ' 



(naive) . 



This is incorrect, for a variety of reasons. In the remain- 
der of this article we will delineate the correct procedure 
to obtain the effective description. 



Fluctuations 

Let us now assume that the spatially-homogeneous 
component 4>o(t) satisfies its own equations of motion to 
lowest order in gi and consider fluctuations around this 
background [70] . 

p(t,x) = 0(t,x)-0o(t). 

Although our observable of interest is the empower spec- 
trum which does not usually require gauge-fixing, we 
wish to incorporate the effects of interactions for which 
the inflaton field and metric fluctuations mix together, 
requiring the construction of gauge-invariant quantities. 



We therefore also consider tensor fluctuations. Utilizing 
the ADM formalism, we parametrize the metric as 

ds 2 = -N 2 dt 2 + hij{dx l + N i dt){dx j + N ] dt) 

with TV the lapse function and N l the shift vector. Sub- 
stituting this into the action produces 



S 



dtd 3 WhN 



M, 



pi 



R 
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2N 2 

2N 2 
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h^di&djcf) (5) 
h^dixdjX-V^x) 



where 



h 
R 



det hij , 

Ricci curvature of spatial metric, 
1 



e vj = ; 2 [h,, Y,X; 



E 



Ei 



The equations of motion for N and N l are just the Hamil- 
tonian and momenta constraints, 



2 
1 



R 



2N 2 



(Eij&t - E 2 ) 



- N% 
2N 2 



(6) 



Cv — N l dv) 2 1 



= 



E 3 l -6/E 
N 



(7) 



N 



N 



In solving for N and N l we can ignore the effect of x; as 
we assume a single-field inflation model, rather the the 
hybrid-inflation model of <p and X- This is valid because 
the corrections to the free equations of motion will be 
0(g 2 ). Since this is the same order as the fluctuations 
we are computing, we may trust our unperturbed back- 
ground solution. 

There are two choices of gauge we will employ: spa- 
tially flat and uniform density. That is, using different 
gauges we may exchange perturbations in <j> for (certain) 
perturbations in g^ and vice-versa. Both will be useful 
at different points in our calculation. 



Spatially Flat Gauge 

We will perform the interaction calculations inside the 
horizon using the spatially flat gauge. This corresponds 
to the choice 



"ij 



a(t) 2 hij, 
0, 
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Fluctuations are then parametrized by ip, 7^ and x- We 
now need the solution for N, N l in terms of these fluctua- 
tions. Since we are interested in quadratic fluctuations of 
ip and 7,-j, it suffices to compute the background correc- 
tions to first order. This is because any third order terms 
for the background would multiply the zero-order terms 
for the fields, which are automatically satisfied. We de- 
note N = 1 + Ni so that Ni is a first-order perturbation, 
as is N i . The solutions are given by 



N 



1 



ei tp 



2 M pl 

Ni = diii, did l n = 



ei d ( 1 <p 

Substituting this back into the action ^ yields (to lowest 
relevant order) 



S 



dtd^x a fi 
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a 3 dt 



(a*Hei] 



1 



{d X f 



1 



y¥> X + yVX ' 



Note the appearance of a small inflaton mass induced 
from the gauge-fixing. In writing this we have omitted 
the ipx mass-mixing term which is subleading in H/M. 
Additionally, we have neglected all interactions which do 
not contain Xi since we are only interested in the ef- 
fect of this heavy field. One of these neglected terms, 
the h 1 ^ diipdjip coupling, produces an M p i-scale interac- 
tion and it would be interesting to compare this to the 
result obtained here. We save this for a future study. 



Uniform Density Gauge 

At the moment of scalar fluctuation horizon crossing 
we then convert these to the uniform density gauge since 
these will then stay constant and are good observables. 
This gauge corresponds to the choice 

hij = a{tf [(1 + 2C)<% + 7y] , 



9 



0, 



daij = 0, 



7r. 



0. 



Note that the x an( i 7 fluctuations remain identical, while 
scalar fluctuations are now parametrized by £. The con- 
version between gauges is 



c 



(8) 



ei M p i 

Of course, were we to include the effects of x on the back- 
ground evolution, this hybrid-inflation model would not 
have constant scalar fluctuations. But as addressed pre- 
viously, the leading-order ^-fluctuations employ a single- 
field background and hence these fluctuations will still 
have constant superhorizon 



THE IN-IN FORMALISM 
In-Out Amplitudes Versus In-In Expectation Values 

Quantum field theory in a static background most of- 
ten employs the "in-out" formalism to produce scatter- 
ing amplitudes. Defining the states |in) and |out) in the 
asymptotic past and future, respectively, amplitudes are 
defined as 



Ac 



(out I Ol in) 



Cross-sections are then obtained by squaring the ampli- 
tude. For non-equilibrium systems, such as a cosmolog- 
ical background, the fundamentally sound approach to 
computing expectation values such as the power spec- 
trum is the Schwinger-Keldysh approach [55 i. 

The procedure is the following. At some early time t; n 
(in the present context, the onset of inflation) we begin 
with a pure state |in(ti„)), then evolve the system for the 
bra- and ket-state separately until some late time t, when 
we evaluate the expectation value: 



(0(t)) = (m(t)\0(t)\m(t)) 



(9) 



Traditionally, the in-state |in) is taken to be the Bunch- 
Davies vacuum state [55], but this is not necessarily so. 
Expanding cosmological backgrounds allow for a more 
general class of vacua, which can be heuristically consid- 
ered to be excited states of inflaton fluctuations. In the 
present context, we will find that integrating out high- 
energy physics generically results in boundary terms in 
the effective action, which represent such excited states. 

If we denote the fields representing the "evolving" ket 
to be {if+, X+i an d those for the "devolving" bra 

to be {<£>_, X-i7ij,-}i the in-in expectation value (JoJ) can 
be computed from the action 



S = S[(p+,x+,lij,+] ~ S[<p-,X-,lij-] 



(10) 



together with the constraint that <p+(t) = ip-(t), X+W — 
X-(t) and 7y,+ (t) = Jij.-(t)- It is then helpful to trans- 
form into the Keldysh basis, 



V = 0+ + ¥>-)/2, $=(yS + -(^_, 

X = (x++X-)A X = x+~X-i 

% = (7«,+ + 7y-)A r = 7lJ - + - 7lJ - 



5 



In this basis the total action ( 10 1 equals 

s[£,$,x,x,7,r] = 

- J d 3 xdt a(t f [-MffiijdFij + d(pd$> + j$ 

2 d 
a?dt 

, 9i 



T 



n $ 2> 
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X 2 
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<^^^^^^^^^^^^^^ ooooooooooo 



1, 



jijdixdjX -r -x „ 



(11) 



Density Matrices 

Although we started with a pure state in ([9]), in general 
we could take expectation values with respect to a mixed 
state, 

(0(i)> = J^imMO^lmjit)}^. 

Here the density matrix p is normalized so that Trp = 
J2iPi,i — 1- l n the path integral language, the density 
matrix is equal to the logarithm of the imaginary com- 
ponent of the action, so that 

S = Rc6> — i ln p. 

We will find that integrating out x wn l generically lead 
to a mixed state for ip and 7, since we are removing states 
from a unitary process. 
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FIG. 1: Feynman rules for our theory in the in-in formalism. 
Single solid lines indicate contractions of <p, dashed single lines 
indicate those of <&, with analogous notation for double lines 
indicating the heavy field components {x,X}, wiggly lines 
indicating the graviton component 7 and circled dots for F. 



Perturbative Solution of Fluctuations 



r defined by 



The equations of motion for the fluctuations are 



8?(p + 3H0-j~[m 2 --^- (a 3 He,) ) <p = 9l p X + ^*X, 



d z X + 3H X - M Z X = v lip 



2 _ _ 9i ( -2 , * 2 \ , 92 



-$X 



4/4 



-9j [ HjdiX+ ^TijdiX ) , 



if 2 - , , :iH% = ^f- Uxdjx + \diXdjX j (12; 



We will now solve these perturbatively. The Feynman 
rules are summarized by Figure [T] where we have used 
the same diagrammatic notation as in [5 7) . 



ds 2 = a(r) 2 (-dT 2 + dx 2 ). 



Neglecting the interactions in equations (12 1, then 



Fourier transforming into the comoving momentum ba- 
sis, they become 



2 d 



^0)» +2H _(0),_ (fc 2 +M 2 a2) -(0) = Qj 



7, 



(0)// 



2.-,(0) 
«7,k 



0. 



Here we have introduced the conformal Hubble parame- 
ter H(t) = a! j a 2 and suppressed all indices on 7. The 
solution for Cp^ is given by the Hankel function of the 
first kind, 



Zeroth Order 

We first consider the free fluctuation equations of mo- 
tion. To do this it is helpful to switch to conformal time 



Za(T) 

~ 3 ^ ^ t2 
2 2 



(13) 
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If we expand in terms of the slow-roll parameters, this is 
the familiar expression 



Mr) 



H 



(1 - ik T )e~ lkT + 



The linearly independent solution U£(t) is simply the 

(2) 

complex conjugate and hence uses He, . 

For the massive fluctuations the free field solution is 
also a Hankel function but can written more transpar- 
ently using the WKB approximation as 



Vk(r) 



1 exp 



M 2 
H(t') 2 -, 



o(r) 



M 2 
H(t) 2 t 2 



1/4 



(14) 



The time-dependent 
U (r) = 



k 2 + 



y/2(k 2 - 

frequency here is 
H ^fy2 T 2 , and the WKB approxi- 



mation |w|/w 2 <C 1 is always valid for H/M <C 1. 

For ^-fluctuations the free solution is again given by 
Hankel function but of a different order. 



2a(T) U " { h 




FIG. 2: Inflaton fluctuation corrections (p' 1 * induced by in- 
teractions. Crosses indicate the zeroth-order solution (p^ 
or x (0) - 




FIG. 3: Statistical averaging of a tree-level solution may pro- 
duce loops, despite remaining classical. Note that the second 
diagram can always be cancelled with an appropriate coun- 
terterm. 



Higher Order 

We may then iteratively solve for higher-order solu- 
tions to the fluctuations, 

V {1 \x) = J d* yy /gG R (x,y) (W°¥ 0) 0/) + f X ( %) ; 

(15) 

Fourier transforming into comoving momentum, these 
vertices can be evaluated using the retarded Green's func- 
tion G R can be written in terms of the fluctuation solu- 
tions, 



G R (n,r 2 ) 



i<4 0) (n)$^(r 2 )) (16) 
-20(n - 75)Im [f/k(n)^(r 2 )] . 



These are shown in Figure[2j The advanced Green's func- 
tion G A is then simply the time-reversal of this: 



G^^G^n). 

Note that comoving momentum is conserved at ver- 
tices. A similar procedure applies for yP^ and using 
their corresponding retarded Green's function Q R and 



Sk(n,r 2 ) 


- ^(xL 0) (n)x^(r 2 )) 


(17) 




= -29(T! - T 2 )Im [V k (n)V£(T 2 )} , 




G£(ti,t 2 ) 








- i<7 k 0) (n)rL° k (r 2 )) 


(18) 




= -20(ri-T 2 )Im[Wk(ri)Wi(r 2 )] 




^(ri,r 2 ) 







Of course the solution may then be iterated again to 
obtain yet even higher-order solutions, 



^ 2 \x) = / dS^g G R (x,y) (W^W (y) + 



Note that for fluctuations the slow-roll parameters are 
implicitly included, and that the order of the solution 
refers to the order of the couplings gi. 



Statistical Correlations 

After obtaining (p to our desired order, we may then 
take statistical averages of the zeroth order solutions via 



F k (n,r 2 ) 





^k(n,r 2 ) 





Mn,r 2 ) 



<^ 0) (n)^ k (T 2 )) 



:(0), 



= 



ri)x^ k (r 2 )) 
Re^n)^)], 

(7 k 0) (n)7i° k (r 2 )> 

Re[W k (n)W£(T 2 )}, 
(rf(n)rL° k (r 2 )) 



to get correlations. This can be heuristically thought 
of as gluing the crosses together in all possible ways. 
An example is shown in Figure [3j While this may pro- 
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FIG. 4: Correction to the background linear coupling coming 
from x-Auctuation backreaction. This can be absorbed into 
the definition of j. 



POWER SPECTRUM EVALUATION 

Here we will give the explicit formula for computing 
the scalar and tensor power spectrum in our high-energy 
model. 



Inflaton Fluctuations 



duce Feynman diagrams with loops, representing inte- 
grals over comoving momentum, they are nonetheless 
completely classical. An important difference of the in-in 
formalism is that loops may represent statistical, but not 
quantum, fluctuations. 



Self-Consistency of Background Solution 

The introduction of fluctuations means that the back- 
ground solution will be slightly modified due to backre- 
action. This is an 0(52) correction to the linear coupling 
j, as shown in Figure |4j having the value 



rf 3 q 92 

(2vr) 3 2 



y (0) y 

Aq A 



(T)> = 



52 
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d 3 q 

(2tt) 3 
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The apparent time-dependence is illusory, as follows. 
Let us convert the integral over comoving momentum 
into one over physical momentum. In static-background 
QFT loop regularization one imposes a UV cutoff A on 
the Wick-rotated 4-momentum, thus respecting Lorentz 
symmetry. Since the quasi-de Sitter background of infla- 
tion breaks Lorentz symmetry ('energy' is ill-defined) we 
cannot impose such a cutoff, but comoving momentum 
is conserved, and therefore we can perform an equivalent 
procedure as follows. Via (14) the X" ne ld's conjugate 
variable to time is given by 



E 



y/q 2 r 2 H 2 



M 2 . 



Thus we may place a cutoff on E at some scale, effectively 
imposing a cutoff on q. Let us momentarily assume that 
H is constant, allowing us to estimate the correction as 



.92 /VA^MV^M d 3 q H 2 T , 

2 Jo 
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M 2 

H 2 T 2 

M 



A + VA 2 - M 2 



We will need to know the correlation of inflaton fluc- 
tuations at the moment of horizon crossing, 

fc 3 

To evaluate this perturb at ively, substitute the classical 



solution (p — (p 1 - ' + (p^ > + ■ ■ ■ into this expression. In 
the decoupling limit Qi —> or M — > 00, the inflaton 
fluctuation power spectrum is simply 



4 0) (fc) 



fc 3 

7T 
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. 2M M\T k )V [ X{T k )) 



F k (T k ,T k ) 



where r k is the horizon-crossing time of mode fc, 
r k ~ — 1/fc. The first order corrections are then: 



There are no interactions in the action (111 which will 



produce this, and so we turn to the second order contri- 
butions, 

P^k) = ^((^l)+2(^)). (19) 

There are a total of four diagrams which could possibly 
contribute at O(gf), shown in Figure [5] We will give 
explicit expressions for the first two, 



d 3 q 



2ir 2 



;' Tfc f Tk 

{-Wif / dri a(n) 4 / dr 2 a(r 2 ) 4 x 



(2tt) 3 



i-«T t ,Ti)] 7- q+k (ri,72 ! ) J F q (7i,7a) [-iG£(T 2 ,T k )] , 



P^ ] (k) = ^ (- tgi f £ dn a(n) 4 £ dr 2 a(r 2 ) 4 x 
/ (2^ H G k(r fe ,n)] H^(n,T 2 )] [-#£(71,7*)] [-<(r 2 ,r fc )]. 



This is independent of H, assuring us that it is the same 
answer one would get from a static spacetime answer had 
we simply truncated E < A. This time-independent cor- 
rection can then be cancelled by appropriate redefinition 
of j. A similar procedure can of course be done for the 
tadpole of x, as shown in Figure 3. 



Inflaton Vertex Evaluation 

We first dissect the interaction at the vertex before 
evaluating the full diagram. Since the time coordinate is 
integrated over, vertex evaluation is very different than in 
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FIG. 5: Power spectrum corrections mediated by the heavy 
field. Single solid lines indicate contractions of if>, dashed 
single lines indicate those of <&, with analogous notation for 
double lines indicating the heavy field components {x,X}. 



standard quantum field theory. Writing out the Green's 
and Wightman's functions in terms of C/'s and V's, we see 
there are three types of cubic vertices. We will evaluate 
each. 

The first is 

A, (r k ; ki , k 2 ) = £ dr a (r) 4 U kl (r ) U ka (r) V* (kj +ka ) (r) . 

The superscript k indicates there is some dependence 
upon the moment that the mode leaves the horizon. Since 
we assume all interactions happen well inside the horizon, 
a good approximation is 



1 



a(r)V2A; 



exp 



-i I kr + 



k\r\ » 1 



where v is given by ( 13 ). To leading order in the slow-roll 



parameters, the scale factor is given by 

1 + ei 



a(r) 



tH 



Substituting this into ( 20 1 results in the expression 

(1 + ei)e-** 



A(Tfe;ki,k 2 ) 



r/r 



where r* is some reference time and -H* is some the Hub- 
ble scale at this time. Then 



du 



H 

M 



(l + ei)dr. 



In this variable it is clear that A\ admits a stationary 
phase approximation at the energy-conservation moment 







d<j> _ M/H 



du 



l + eidr 



M/H 



du , , 
7T— + ki + k 2 
dr 



1 + ei 

The variation in slow-roll phase is given by 
dv dN dv 3 



ki+k 2 | 2 



dr dr dN 



and so is higher-order in slow-roll and can be neglected. 
The solution then defines the NPH. 



- 1 = - V / 2M^(I - 



cos( 



cos( 



ki ■ k 2 

k\k 2 



Expanding the phase near this stationary point gives the 
inflection 



d 2 ^ 
du 2 



M/H 



1 + ei 



ki + k 2 



Then to leading order in H/M and slow-roll parameters 
the amplitude is 



-4i(T fe ;k 1 ,k 2 ) 



(22) 



Triflfa - t c )9(t c - r in )(l + |ei)e-^ 



2y/hk2 [2fcifc 2 (l - cos(9)] 1/4 VHM 



2V2hk~ 2 J Tm Hr {]ki+k2]2+ m^ } i/4 

(21) 

where the phase is given by 



The physics of this is clear. This diagram accounts for 
the threshold production/decay of heavy particles at high 
rcdshift in the early universe. Although this induced 
term is localized in time, it is fc-dependent and thus does 
not represent a physical boundary. Rather, it represents 
a 'hypersurface' representing the energy scale of interac- 
tion, precisely as expected in the New Physics Hyper- 
surace approach. 

The second possible vertex is identical to Ai but with 
one U conjugated, 



(r) = ^ + (fc 1 +fc 2 )r- / dr'W|k 1 +k 2 | 2 + -— -. p 

Jn n V H 2 t u A 2 {T k ;k 1: k 2 ) = J dr a(T) 4 f/ kl (r)C/^ 2 (r)V"_! (kl 



,(r). 



Let us now introduce the rescaled time coordinate u such 
that 

H 

M 

To obtain the differential, recall that the Hubble scale 
changes with time as 



H(t) ps H* 



This has only imaginary-time saddlepoint solutions. 
Since our r-integral is confined to the real axis we will 
never pass over this point in our integration, and so 
this amplitude will be suppressed as A 2 ~ erf(^) ~ 
H_ e -(M/H) exce p^ f or cos « —1, when the saddlepoint 
moves to infinity and the inflection vanishes. In all other 
cases we may neglect such interactions. Finally we con- 
sider A3 which has both t/'s conjugated and so admits 
no saddlepoint solutions, and thus can also be neglected. 



9 



1-cos 



tfin 



FIG. 6: The phase space of diagram A in both {q, 1 — cosd} 
and {r, E} coordinate systems. Constant-r contours arc 
shown in blue, and a constant- E contour is shown in red. The 
low-energy phase space allowed in our integration is shown in 
gray. 



Inflaton Power Spectrum 

After the vertex evaluation it remains to perform the 
integration over loop momentum. 

First scalar correction 
Let us first consider diagram A, 

,2 j.3 _ r ^3 



(27T) 



C1 ^(r^q)! 2 . 



. (23) 

This is non-trivial for two reasons. First, we wish to 
implement a physical cutoff, and so a change of coordi- 
nates is required. Second, loop momentum influences the 
time of vertex interaction via r c (q, k), and so affects the 
Hubble scale. Thus we need to convert the integral over 
comoving momentum to one over time of interaction. To 
accomplish this, first simplify the 3d integral to 2d, 

Now replace the coordinates {q, 1 — cos 9} with {t, E} via 
the definitions 



H 

M 



E = Hq\r\ = M 



y / 2kq(l - cos9), 



2fc(l 



Of course the definition of r coincides with r c , since this 
is the stationary-phase localized time of interaction. In- 
verting these coordinate systems gives 



The measure transforms as 



d(l — cos 9)dq — > 



M 2 {l + e 1 ) 
H 2 Er 3 k 



dEdr. 



The integral in ( 23 1 is then transformed as 



' / " q |./4i(7*;k,q)| a 



(2tt) 3 
1 



^/qdqd(l — cost 

16V2nk 3 / 2 M J Hy/l-co&6 

1 f drdE „ , 
(l + 4ei). 



-(l + 3ei) 



16nk 2 



The phase <fi cancels out in this amplitude and so we 
need not evaluate it. We can now easily place limits on 
the region of integration in terms of the energy scale A: 
Maximal Q As done previously, we place a limit on 
the physical energy to be A. Using the fact that the 
vertex with a virtual ip and an external if of comoving 
momentum k have precisely the energy of the %-field, the 
bound is found to be 



EL 



-HkW\ + A. 



Minimal E 
limit on E, 



Since 1 — cos < 2 this places a minimal 



M 2 



AH(T) 3 k\i 



Although this is time-dependent, it remains scale- 
invariant because it contains only the combination kr. 

Initial and Final t Of course we require 
Emm < Emax- Their equality will determine the 
earliest and latest times for which this condition is valid, 
an equation which is quadratic in r. The two solutions 
are given by 



.(0) 



JO)' 



1 - ei In 



where we have defined 



(o) A ± VA 2 - M 2 



2hH* 



(25) 



The later root is 



AO)' 



Tfin =t! 1 - ei In 



whereas the earlier root is 



-| 0; [ 1 - ei In 



(0)' 



This will always be after the point where the external en- 
ergy parameter E^ = Hk\r\ reaches A, so we are assured 
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that all energy scales under consideration are below A. 
Figure [6] shows the phase space in the two coordinate 
systems. 

This gives us a window T; n < r < Tfi n in which to inte- 
grate out this high-energy interaction in a controlled way. 
Note that choosing A < M makes this window vanish, a 
reassuring fact since (by definition) there should not be 
any interactions below the scale of New Physics. There 
is an upper bound on A as well: the interactions must 
cease before the gauge-transformation at the horizon-exit 
near = —1/k (whereupon our plane- wave approxima- 
tion would break down). This means we are allowed to 
set the cutoff scale in the window 



M < A < - (H + M 2 /H) 



For a typical value of M /H ~ 100, this means we can set 
M < A < 50M + H, which should be adequate for the 
low-energy behavior. We now expand the Hubble rate 
around the reference point In r* = — In fc» , and keeping 
only the leading term in A/M, the integral is evaluated 
as 



P. 00 



-3 I 1 + 4ei - 3ei In — 



8A 3 H 2 k 2 
3M 4 



1 - -ei + 4ei In 



dE 
M 2 



The power spectrum correction is then 



4 A) (fc) 



fit 



2 24tt 3 
2 



\U k (r k )f 



1 



A 3 
fAAH„k 



3 £1 4£lln \~APk 



(26) 



Second scalar correction 

Diagram B can be evaluated in an identical fashion, 
except for two subtleties: the Heaviside function in the 
Green's function allow only half of the Gaussian fluctua- 
tions to contribute, and the two i's produce a sign change. 
Thus 



In k 



FIG. 7: Power spectra for the example theory in slow-roll 
inflation. The red line shows the effect of interactions relative 
to that of the flat spectrum of the free theory in exact de Sitter 
space shown in black. Note the corrections slightly decreases 
at higher k as a result of the Hubble scale changing. 



Total result 

Summing diagrams A, B, C and D results in the fol- 
lowing total shift in the scalar power spectrum, 



AP v (k) 



Silk 3 
96tt 3 



A 3 

2 fAAH^k 



(27) 



Since A ~ M and g± ~ H , the correction then scales as 
AP ~ H/M. Figure [7] shows these corrections to the 
power spectrum in the example theory. 



The Tensor Power Spectrum 

We did not include a 7 2 % coupling in our high-energy 
model, nor are there any induced by gauge-fixing; hence 
the tensor spectrum is the same as the standard slow-roll 
case. The gravitational power spectrum is 



p(0) 

r h,ij 



(k) = 



32k 3 
32k 3 



5yF(Tfc,Tfc). 



(28) 



p(B) 



(k) 



Comments and Interpretation 

0, the result is a perfectly 



In the de Sitter limit of e 
scale-invariant correction, 



Third and fourth scalar correction 

Diagrams C and D can also be evaluated in a similar 
fashion, and to leading order in H/M they cancel. 



9P 



192tt 3 M 4 ' 



Let us verify that we should have anticipated this by 
simply computing the physical momentum p at the NPH 
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FIG. 8: The evolution of fluctuations during inflation, as they 
encounter the NPH. The scale of New Physics M is again in- 
dicated, as well as the cutoff scale A. The comoving momen- 
tum of a fluctuation is given by its slope. The high-energy 
interactions in the blue region can be effectively described as 
excited initial states for the inflaton. The fact that all fluc- 
tuations encounter New Physics at (nearly) the same energy 
scale ensures that the answer is (nearly) scale-invariant. 



for some fluctuation of comoving momentum k: 

p(r e ) = k/a(r c ) = -H(r c )kT c = M/2. 

So the NPH is indeed at a fixed energy scale for all modes. 
A second verification is to compute the period of inter- 
action near this NPH. After all, we used a stationary 
phase approximation which includes the gaussian fluctu- 
ations around the stationary phase, giving the NPH a 
finite width. Re-expressed in terms of physical time At 
shows this interaction period also remains scale-invariant: 




Figure [8] shows a schematic of this. 

The dependence upon H/M is also worth comment- 
ing upon. Recall the argument presented by [3H5TJ] in 
considering the low-energy expansions of local operators, 



Sb.uik = / d 4 x 



Fluctuations frozen at p ~ H will then have a power 
spectrum scaling as 



H 2 1 



XH 2 
~M 2 



displaying high-energy corrections at order (H/M) 2 . The 
addition of boundary terms such as 



5, 



boundary 



= d 



X^ 



d l ipdiip 



then modifies the power spectrum at order H/M, poten- 



tially a more significant observable effect. Clearly, (27) 
does not reflect this expectation. Where did the argu- 
ment fail? 

The answer lies in the assumption of a local operator 
expansion. The operators in a UV-complete theory are 
of course local, and this property is also respected by the 
low-energy (p 2 /M 2 <C 1) expansion of massive Green's 
functions in static spacetime: 



G(x-y) = - 



d 4 p e -^-^-y) 
M 2 



1 

M 2 



(2tt) 4 p 2 
1 + 



M 2 



5\x-y). 



The Green's (and Wightman's) function used in quasi-de 
Sitter backgrounds, however, do not yield local effective 
forms. Instead, the appropriate expansion parameter is 
H/M <C 1 which we have seen via the stationary phase 
approximation is an approximation to energy conserva- 
tion. In this expansion, non-local operators are generated 
which betray the expectation of the small momentum ex- 
pansion. This was explicitly shown in |69j . 



Parametrizing Power Spectra 

Using the gauge-conversion in eq. (JsJ) , the scalar per- 
turbation is then 



P d k) = i^. 



(29) 



The advantage of using comoving k to parametrize the 
power spectrum is that there is a well-defined map to the 
CMB ^-modes. To compare such predictions to observa- 



tion requires summarizing the power spectrum ( 29 ) by 



a few parameters which can then be constrained by the 
data. The late-time limit of the modefunction is 



Hm 



l-ex-lei 



C-Ui 



2 / nJ* 



Here we have used the abbreviation 

r H3/2) . _ 

C -rP)-' n2 

= 7 + In 2 - 2 

where 7 is the Euler-Mascheroni constant. Following 
the example of Martin and Ringeval [58], we can then 
parametrize the scalar power spectrum as: 



Pco(h) 
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where 



(30) 



+ 



a 



{ s) = 1 - 2(C + l)e x - Ce 2 
8ei 
3 

5?A 3 



48ttM 4 H« 



1 



-2ei - e 2 



2ei + e 2 )C - 4ei In 



(6ei + e 2 ) . 



/4AiJ» 



24nM 4 H a 



While the search for oscillations in the WMAP3 data [55] 
failed to show evidence for such high-energy oscillations, 
the Planck satellite data will be of much higher resolution 
and is worth investigating. We discuss this matter further 
in the concluding section. 

For the gravity waves, which have no explicit high- 
energy interactions, the spectrum reduces to the standard 
slow-roll result. The late-time limit of the modefunction 



l-ei((7+l)-eiln- 



The parametrization is then 



where 



Ph(k) 
Pho(k*) 

Pho(h) 



16ff? 



a 



(t) 



= l-2(C+l)ei, 
- -2ei. 



OBSERVATIONAL POSSIBILITIES AND 
CONCLUSION 

In summary, we have developed a technique to explic- 
itly calculate the generic corrections to the inflaton power 
spectrum from fundamental high energy physics. These 
can be connected to microscopic models through a "gen- 
eralized" Boundary Effective Field Theory as in [38] . 

Why does the naive effective action of the power spec- 
trum fail? For three reasons: 

1. Correlation functions (such as the power spectrum) 
must be computed using the in-in formalism, sig- 
nificantly altering the role of classical and quantum 
fluctuations. 

2. The expanding background means that heavy fluc- 
tuations can now be excited on-shell. 

3. The interactions producing such heavy modes 
(whether on-shell or virtual) are now localized at 
specific moments in time, giving rise to effective 
interaction boundary terms. 



Because there are no distinct features in the power 
spectrum, any visible new physics will probably need to 
be found in higher order correlations. The bispectrum, or 
3-point correlation [63J , for this example model will van- 
ish to leading order in slow-roll, but more general models 
should have the effects of inflaton self-interactions and 
modified vacua multiplying |64L I65| . Our results justify 
the conjecture made in [65] that effective actions and 
modified vacua are expected in tandem. One distinct 
hope is that the bispectrum may be oscillatory. 

For the trispectrum, or 4-point correlation [68 , lit- 
tle has been done in the way of categorizing such cor- 
relations and we cannot give a complete analysis here, 
but let us summarize a few basic ideas. Our model has 
leading-order trispectrum terms of types ((p 2 )((p 2 ) and 
(</? 4 ). These include the characteristic oscillations which 
are now squared, and would display high-energy signa- 
tures even if the oscillation frequency ~ H/M were too 
small to be observed by the power spectrum [71] . 

Looking towards complete understanding, further de- 
tails of the cutoff procedure, renormalization group prop- 
erties and vacuum selection in an expanding background 
need to be well-understood. 
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